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Abstract — We consider a hybrid wireless sensor network with 
regular and transmit-only sensors. The transmit-only sensors do 
not have receiver circuit, hence are cheaper and less energy con- 
suming, but their transmissions cannot be coordinated. Regular 
sensors, also called cluster-heads, are responsible for receiving 
information from transmit-only sensors and forwarding it to 
sinks. The main goal of such a hybrid network is to reduce the 
cost of deployment while achieving some performance constraints 
(minimum coverage, sensing rate, etc). 

In this paper we are interested in the communication between 
transmit-only sensors and cluster-heads. We develop a detailed 
analytical model of the physical and MAC layer using tools 
from queuing theory and stochastic geometry. (The MAC model, 
that we call Erlang's loss model with interference, might be 
of independent interest as adequate for any non-slotted; i.e., 
unsynchronized, wireless communication channel.) We give an 
explicit formula for the frequency of successful packet reception 
by a cluster-head, given sensors' locations. We further define 
packet admission policies at a cluster-head, and we calculate 
the optimal policies for different performance criteria. Finally 
we show that the proposed hybrid network, using the optimal 
policies, can achieve substantial cost savings as compared to 
conventional architectures. 

I. Introduction 

In this paper we analyze performance of a hybrid sensor 
network architecture, proposed in [1]. The goal of this archi- 
tecture is to reduce the cost of deployment while achieving 
some performance constraints. These constraints concern the 
event-to-sink performance of the network [2], and thus require 
a sufficient density of deployment to correctly monitor the 
sensing domain, and efficient transport solutions in order to 
provide the detected information to the central unit. 

The idea proposed in [1] assure some fraction of sensing 
capabilities by the transmit-only sensors, who do not have 
receiver circuit, hence are less energy consuming, and cheaper 
[Rabaey, ODonnell, Sheng]. The remaining part of the sensing 
and the totality of the information transport task is assured by 
regular sensors, also called here cluster-heads. In particular, 
they are responsible for receiving information from transmit- 
only sensors, who send it blindly, and forwarding it to sinks. 

An immediate consequence of the assumption on the 
transmit-only sensors is that the transmission traffic generated 
by them is completely random, and thus, we have to admit that 
some part of the detected information will be lost because of 
the collisions generated at the receiver of the cluster heads. 
Note that arbitrarily increasing the density of the transmit- 
only sensors and/or their traffic intensity we may saturate the 
transport layer and thus make the situation even worse. 
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In [1] the authors propose a simple heuristic of a packet 
admission policy at a cluster-head in order to maximize the 
total number of captured packet. In this paper, we define 
a detail mathematical model of the hybrid network. Using 
this model we prove that the heuristic from [1] is indeed 
optimal. We also derive the optimal policy that maximizes the 
coverage region of a hybrid network, which was previously 
unknown. Finally, using our model, we are able to quantify 
the substantial savings obtained with the hybrid architecture. 

In addition, the MAC model that we develop and call 
Erlang's loss model with interference, might be of independent 
interest as adequate for any non-slotted; i.e., unsynchronized, 
wireless communication channel. 

Related work: The only work on transmit-only sensor 
networks we are aware of is [1]. An analysis of the event- 
to-sink performance of standard network is given in e.g. [2]. 
Several works show a significantly lower complexity and 
power consumption of transmiter over receiver circuit [3], [4], 
[5]. 

Note that we are not interested in this paper how to obtain 
the sensor deployment that satisfies a correct monitoring of the 
domain. Several results on this subject were already published; 
e.g. using the explicit formula for the volume fraction of 
the stochastic geometry Boolean model (see eg. [6], [7]) or 
asymptotic formula, which can be applied when the density 
of nodes is large while the sensing ranges are small; see [8], 
[9]. 

The remaining part of this paper is organized as follows. 
In Section [11] we describe the system assumptions. Next, 
in Section |lll] we evaluate and optimize its performance. 
Numerical examples are presented in Section Hvl In Appendix 
we develop some details of the mathematical model used to 
analyze the system. 

II. System assumptions 

Let us consider a network of sensors and cluster-heads. 
Sensors are simple sensing devices that are equipped only 
with a single transmitter They are supposed to sense and 
periodically send information to the cluster-heads. Cluster- 
heads are more powerful (and more expensive) sensors. They 
are equipped with a receiver and a transmitter, and their special 
role is to collect information from transmit-only sensors and 
forward it to a central server. The network consists of a large 
number of sensors and a much smaller number of cluster- 
heads. We want to analyze and optimize the performance of 
the information transport from sensors to cluster-heads. 

A. Events and traffic 

We assume that some events trigger transmissions at the 
sensors randomly and independently of each other, with in- 
tensity Ag. Two scenarios are possible. In the first one, event 



is a time instants at which a sensor decides to transmit 
information about the actual state of the sensed environment. 
There, Ag is a system variable controlled by the designer. In 
the second scenario, event is a time instant at which a sensor 
senses some random excitation in its proximity and transmits 
a report on it. A local character of random excitations justifies 
the assumption of independence between transmissions of 
different sensors. In this case, Ae is an external parameter 
depending on type of events we measure. 

Our analysis applies also to a situation when the excitation 
of the medium is not local and persists for some time (like 
in an intrusion detection problem). Then, the time and space 
scale of our spatial throughput analysis corresponds to the 
duration of the excitation and the region from which the 
sensors report about this excitation. A random. Aloha-type, 
back-off mechanism has to be implemented at sensors, in 
order to avoid systematic packet collisions from other sensors 
sensing the same excitation. The vent is a time instant at which 
the back-off mechanism of a sensor makes it transmit; thus, we 
find the first scenario in this interpretation. When the excitation 
is over, the sensors may go to a sleep mode, or sense at some 
smaller rate. 

A canonical example, considered in this paper, is this of a 
channel with the throughput of 1 MB/s at the physical layer, 
which is shared by sensors homogeneously distributed with 
the density of As = 10 sensors/m^. We consider the periods 
when the channel is actively used by all the sensors which emit 
with the temporal intensity Ag = 1 kbps. Such a relatively high 
intensity may be reasonable to perform correctly during some 
periods of persistent excitations. 

B. Reception 

Since sensors are transmit-only, they cannot sense collisions, 
and send packets blindly. The goal of cluster-heads is to 
receive these packets. We suppose that a packet is correctly 
received if the SINR, empirically averaged over the reception 
duration, is higher than some threshold. Otherwise, the packet 
is lost. We consider a slow fading Gaussian channel channel 
model with repetition coding with interleaving. Repetition 
coding corresponds to CDMA or UWB spreading. Interleaving 
means, that each bit is sent through many symbols uniformly 
distributed over the duration of the packet size. Suppose 
that the receive is equipped with a matched filter (coherent 
maximal ratio combiner). Then, the standard analysis (see 
e.g. [10, Section 3.2.1]) says, that there is a threshold 7 on the 
SINR that should be respected in order to maintain the link 
quality; 

'^''^Tm . >7, (2.1) 
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where W and Ij is, respectively, the noise power and the 
power received from interferers during the j symbol, P^ec 
is the received power averaged over fading effects that is 
supposed to depend only on the emitted power and the distance 
between the emitter and the receiver. Remark, that interleaving 
allows for the empirical averaging of the interference / in il.H 
over the packet duration. We will interpret i2.1\ as the SINR 
condition identifying the successful reception of the packet at 



Fig. 1. Gain due to admission policy. 

the MAC layer, given the link fading value h and the received 
interference power process {Ij}- 

Commonly used fading model is Rayleigh fading where 
link fading h is represented with a circular complex Gaussian 
random variable. In this case can be seen as a realization 
of some exponential random variable (see e.g. [10, p. 50 and 
501]). 

Cluster-head communication: We also assume that 
cluster-heads have a reliable communication of a higher rate 
to a central server, which does not interfere with the sensor 
channel (e.g. can be wired, but not necessary). 

C. Synchronization and decoding 

A cluster-head needs to synchronize to a packet sent by a 
sensor before receiving it. In order to synchronize to a packet, 
the cluster-head has to receive it with some minimum reception 
power ^. If the packet reception power is higher than this, 
cluster-head starts receiving and it continues receiving until 
the end of the transmission. // the transmission is lost because 
of the collision with another packet emission (interference), 
the error will be detected only at the end of the reception. 
Moreover, this interfering packet will be lost as well since the 
cluster head was not idle at its arrival epoch; cf. Figure [0 

In order to improve the efficiency, we introduce a packet 
admission policy. Once the cluster-head is synchronized to a 
packet, it may decide to receive or ignore the detected packet 
according to some packet admission policy based, for example, 
on the value of the received power This policy allows to ignore 
some weak packets so that the cluster-head is more often 
available for stronger packets. The choice of a particular policy 
depends on system design goals, described in Section III-EI 



D. Sensor placement 

Sensors in a given network realization are fixed. We show 
how to evaluate the performance of a given fixed configuration 
of sensors served by one cluster head applying some admission 
policy. For performance optimization, we adopt a stochastic 
approach. Namely, we optimize the system parameters with 
respect to the average performance, taken over all possible 
spatio-temporal configurations of packet emissions, which are 
driven by some spatio-temporal Poisson point process. We 
call this model a Poisson-rain process of events. We argue 
in Appendix lA-B. II that this model — which corresponds to 
the situation where each event is associated with one sensor 
that activates only once, when the event occurs, and leaves 
the system afterward — is a reasonable approximation of the 
packet traffic generated be an arbitrary (also deterministic) 
repartition of sensors on the plane, provided the density is 
large. However, it is important to keep in mind that this 
model is just an approximation, that in reality the sensors are 
persistent, and that we can collect and keep some information 
about each one. 

^e.g. needed to have symbol-error rate during the preamble of 10% 
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E. Design goals and performance metrics 

We consider several design goals, related to the transport of 
the information from sensors to cluster-heads. Our principal 
performance metric is the spatio-temporal density p{x) of the 
received information. We define it as the mean number of 
packets received from sensors, per second, from the surface 
area dx. 

In order to obtain some desired density p( ) the system 
designer can influence the placement of nodes (at least the 
density), transmission power, and frequency of transmissions. 
However, while the density of nodes can depend on the 
location, the transmission power and frequency is the same 
for all the sensors as we suppose we cannot reconfigure each 
sensor separately once it is placed. The main parameter that 
can be tuned in order to shape the information transport from 
sensors to the cluster-heads is the packet admission policy 
applied by cluster heads. 

In general, the network designer may be interested in maxi- 
mizing the coverage of some sensing domain or in increasing 
the total throughput. These contradictory goals are realized by 
the policies respectively called max-min and globally optimal, 
which are defined in Sections IlII-BI and IIII-Cl 

Finally, the goal of this paper is to propose a hybrid network 
that will minimize the cost of the network deployment. We 
thus want to show how much money can a designer save by 
combining wireless tranceiver sensors with the cheap transmit- 
only ones. To that respect we study the economic optimization 
problem in Section UlI-DI which finds the right proportion of 
the two types of devices that minimizes the cost of a network, 
under the constraint of some desired level transport-aware 
effective coverage. 

III. Analysis 

In this section we will analyze and optimize the perfor- 
mance of the sensors-to-cluster-heads information transport in 
transmit only sensor network, whose system assumptions are 
described in Section HI] We assume that the density of sensors 
\s is large enough and the sensing regions are small so as 
to guarantee a good sensing-coverage of the domain with a 
sufficient resolution. 

A. Density of information 

First, we study the density of information received by a 
single cluster head p{x). Remind that we define it as as the 
mean number of packets received by the cluster-head per unit 
of time from the area da;. We will model the traffic described 
in Section Hi- Al bv a spatio-temporal Poisson point process of 
events"* (packet transmissions) with intensity As(da:) x Aedi, 
where Ks{dx) = \s{x)dx and Ae are, respectively, the mean 
number of sensors placed at dx, and the temporal intensity 
of packet traffic sent by each sensor. This Poisson rain of 
events (packets) is supposed to be received by one cluster- 
head, whose behavior is described in Sections III-Btill-Cl 

"^We explain in Section IA-B . II that tliis model, called there the Poisson rain 
of events, is a reasonable approximation of the process of packets transmitted 
be an arbitrary pattern of sensors (not necessarily Poisson) densely distributed 
on the plane. 




Fig. 2. Three terms of interference Ci, C2, Cj^. 

Specifically, the cluster head applies some admission policy 
d{x), which is the probability that it tries, given it is idle, 
to receive a given packet emitted from dx ^ . We assume that 
the admission decisions are taken independently of each other 
and of anything else, and thus the spatio-temporal process 
of admissible packets is the Poisson process with intensity 
d{x)\s{x)Ax X Aedi. 

Inspired be the channel description of Section Hl-BI we as- 
sume that a given admissible packet, arriving when the cluster- 
head is idle, is correctly received if some SINR, empirically 
averaged over the reception period B, is higher than some 
threshold 7; cf. ( 12. 1> . The interference is created by all other 
emissions taking place at this time period and by some external 
noise W . A detailed mathematical analysis of the performance 
of the cluster-head modeled by some Erlang 's loss system with 
interference and SINR condition ( lA.ll i is done in Section lA-BI 
under the assumption of Rayleigh fading. In what follows 
we summarize the results of this analysis. First, we remind 
a general fact that follows from the Campbell formula. 

Proposition 3.1: The density of received information is 

equal to p{x) = XeXsix)d{x)pfreePrecix), wherC Pfree IS 

the probability that a typical admissible packet finds the 
cluster head idle when it arrives and Prec{x) is the conditional 
probability that the typical admissible packet arriving from dx 
can be correctly received, given the cluster head stars receiving 
it. 

Suppose that the cluster-head is located at the origin. Denote 
by P is the emission power used be all sensors, by L{x) the 
power attenuation function (path-loss) of the distance from x 
to 0, and hy Cw the Laplace transform of the power W of 
the external (white) noise; CwiO — if this power is 

constant. 

For a given admissible packet received by the cluster-head, 
let £1, £2, C.Jb '^he Laplace transforms of the interference 
averaged over the reception period, generated respectively, 
by: admissible packets arriving when it is being received, 
admissible packets that are being sent at its arrival epoch, 
all non-admissible packets; cf. Figure |2l These Laplace trans- 
forms are explicitly given by formulas ( lA.lU . ( IA.12t . ( IA.13t 
with A = Ae / d(r)\s[x)dx begin the total intensity of the 
admissible packets (the integral is taken over the whole domain 
of the network deployment). Denote — 1 / {PL{x)). By 
Corollarv lA.6l we have the following result. 

Proposition 3.2: The Erlang acceptance probability is equal 
to pfree = 1/(1 + AS) and the conditional reception proba- 
bility is equal to Prec{x) = Cw{lx)^l{lx)^2i-fx)^JB{lx)- 

Lets denote 

£(0 =exp^A,sy(l-^pl^log(l + ePL(x))A,(dx))^ . 

(3.1) 

^The cluster-head does not need to know the location of the receiver; it can 
apply some admission policy depending on the received power. 
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Corollary lA. 71 gives two more explicit bounds on prec{x)- 
Proposition 3.3: We have p^^^{x) < Prec{x) < p^^^{x), 

where p^^^{x) = Cw{lx)C?{lx), Pred^) = Cw{lx)C{lx) 
and £ is given by ( 13. 1> . 

Denote by respectively, the lower and the upper 

bound of p{x) obtained when Prec{x) in Proposition 13.11 is 
replaced by, respectively, p^^^{x) and p^^^(x). Note that the 
both Prgc(^) ^"^d P (^) depend on d{-) which makes 

the analysis of p{x),'p(x) easier (for the quality of the bounds 
see Figure |3j (left)). 

Before describing some optimal policies, we define a naive 
policy dnaive{x) — l{x E Dq). whcrc the set Dq is fixed such 
that mean received power is large enough to receive correctly 
the packet, given only external noise W (no interference) ^; 
i.e.. Do = {a; : PL{x)/W > 7}. 

B. Optimizing the transport-aware coverage 

Knowing that the attenuation function L{x) (and thus 
Prec{x)) typically decreases with the distance to the cluster- 
head, one has to compensate it with an increasing density of 
sensors As(-) and/or a spatial admission policy d{x). 

In this paper we suppose now that the sensors are already 
deployed ^ with some given with density As {x) > on some 
sensing domain D. We look for an admission policy d{x), 
such that any increase of the ratio p{x)/D{x) on some set 
dx of positive Lebesgue 's measure would be at the expense of 
decreasing of some already smaller ratio p{y)/D(jj) on some 
non-null set dy G D. The policy dmaxm{x, D) realizing the 
above principle is called weighted max-min fair policy, with 
weights 1/D{x). It is known that if dmaxm{-, D) exists then it 
is unique. For brevity we will denote by dmaxm (x) the max- 
min policy with equal weights {D{x) = D) (the policy does 
not depend on the value of D). 

We cannot exactly characterize the max-min fair policy 
for p{x), however, we can do this for some bounds. Denote 
/ = Jj^D{x)/p {x)dx. Assume that the sensing domain 
D is compact, D{x) continuous on B and denote M = 
maxxeo D{x) / {Xs{x)p J. ^^{x)). We define in the similar man- 
ner /, M replacing p^^^{x) in the above formulas by p^g^(x). 

For a given poHcy (i(-) denote by ||d||A^ = Jjj,d{x)Xs{x) dx 
the total spatial intensity of admissible packets under d{-). 

Proposition 3.4: • The max-min fair policy 
imaxmi'TD) for p[x) on D exists if and only if 
M < 00, it is equal to 

D{x) 



—maxin 



{x,D) 



(3.2) 



MXs{x)p^Jx) 

and reahzes p {x, D) = D{x) /{Bl+M_/Xe). More- 
over, under L>) we have p{x) > P„^^^Jx, D). 
• The max-min fair policy dmaxm{-,D) for p{x)/D{x) on 
D exist if and only if M < 00, it is given by i3.2\ with 
M. v replaced, respectively, by M ,p^^^, and realizes 
'Pniaxmi^i D) = D (x) / {B7 + JI / Xe) ■ Morcovcr, there is 
no policy d(-) for which p{x) > p„Q^„j(a;, D), with the 
strict inequality on some non-null set dx. 

^This may correspond to the successful synchronization to the packet 
'leving the the optimal deployment problem for future work 



Proof: We consider the lower bound. The proof for the 
upper bound is analogous. Suppose that M < 00. Note that 
the function given by the right-hand-side of J3.2t in positive 
and not larger than 1. Thus d^^^^(-,D) is a policy. Note 
also that for xq — arg max^^j^ D{x) / {Xs{x)p (x)) we have 
dimaxmi^o, D) — 1. Assume now that for some policy d'{x) 
the respective ratio p'{x) > D{x)l/{BI + M/Xe) and that the 
inequality is strict on some non-null set dx. It easy to show 
that then Uri'lU^ > ll^maKmlU., and thus ^'(xo) < p{xo). This 
shows that d^naxm max-min fair 

Suppose now that M = oo. Take any policy Note that 
p{x)/D{x) cannot be constant under this policy (there is no 
such policy). Thus, there exist xi,X2 such that p{xi) > p{x2). 
Note that we can slightly increase d{xi) and decrease d{x2) 
is such a manner that IIcZIIa^ remains constant. This increases 
p{x2) without changing p{x) for x ^ xi, X2. Thus, d{-) is not 
a max-min fair The remaining part of the result follows from 
Proposition 13.31 ■ 

Remark: Suppose the cluster-head is to collect information 
sent by sensors in a given compact domain D with some min- 
imal density p{x) > D{x). The problem might be infeasible. 
However, if it is, policy dmaxmix, D) satisfies the constraint. 

Example 3.5: Consider a uniform coverage D{x) = D x 
l{x E B) witght function. We might be interested in maxi- 
mizing the constant density D given the domain D. This is 
achieved using dmaxm- Alternatively, we might be interested 
in maximizing the area of domain D while providing some 
minimal density D. For example, for a homogeneous repar- 
tition of sensor Xs{x) = Xg and distance-dependent path-loss 
L{x) — L{\x\) model, we maximize the radius R of the disk 
D = -6(0, R) centered at 0, under the contranit p{x) > D 
for X E a. Using Proposition 13 .41 one can find the solution 
R = Emaxm such that poficy d^axm on D = S(0,i?„„^„) 
satisfies p{x) > p{x) = D for all x E O. We illustrate this 
problem in Section Hvl 



C. Optimizing the total throughput 

Consider now the problem of the maximization of the 
total weighted intensity of received information U = 

^j^p{x)/D{x)dx ^ where D{x) > are some arbitrary 
weights. 

Denote by C/, U, respectively, the lower and the upper bound 
of the total weighted intensity of information obtained when 
p{x) is replaced by p{x) and 'p{x). As previously, we can 
solve this global optimization problem for the bounds p{x) 
and p(a;), and in this way approximate the solution of the 
original problem. 

Denote the following water-filling region D(6') = {x S D : 
Xs{x)p {x)/D{x) > 9} and the constant 



arg max ■ 



Io(e) >'s{x)p^^^{x)/D{x) dx 

•^e-B /jj^g^ Xs{x) dx 



m 
1 



We define in the similar manner D, 9 replacing p^^^x) in the 
above formulas by Pj.^^{x). 
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Proposition 3.6: • The poHcy d* [x) — l{x 
maximizes U_. Under this poHcy 

Ae 4(9.) Xsix)p^Jx)/D{x) dx 



Moreover, under d*{-) we have U > U_*. 



(3.3) 



• The poHcy d (x) ^ l{x £ D{9 ) maximizes U_. Under 
this poHcy U — U it is given by J3.3t with D,9_* 
replaced, respectively, by D, 6* , Moreover, there is no 
policy (•) under which U > U . 

Proof: We consider the lower bound problem 
(proof for the upper bound is analogous): maximize 
Jp Xsd{x)p^^^{x) / D{x) dx/{l — \eBA) under the constraints: 
A = Jjj Xs{x)d{x) dx, < d{x) < 0. We write the Lagrangian 



L{d,0, fiQ, pi) = 9A + / fj,i{x)dx 



d{x) 



Xs{x)p(x)/D{x) 



- 6 + pq{x) — pi (x)^ dx . 



1 + X^BA 

By the strong duality and the KKT conditions the optimal 
policy has the form of the indicator function d*{x) = 
l{Xs{x)p^^^[x)/{D{x){l + XeBA*)) <9*) for some 9*, A*. 
The values of these constants are found by the standard watter- 
filling policy. The remaining part of the result follows from 
Proposition 13.31 ■ 
Example 3.7: Consider equal weights D{x) = D, homoge- 
neous repartition of sensors and and distance-dependent path- 
loss model. Then D(0*) = -0(0, i?*) is a dies of radius R* 
We illustrate this finding numerically in Section IIVI 

D. Optimizing the network cost 

Suppose that one transmit-only sensor costs Cs, while a 
transport-reliable sensor (with the same sensing functional- 
ity) costs Cc- Consider an architecture where the transport- 
reliable sensors act as cluster-heads considered previously in 
this paper; call them cluster-heads. Assume that information 
obtained (sensed) directly by cluster-heads sensors is delivered 
to the central unit with probability one, while the information 
obtained by a transmit-only sensor located at x is delivered 
there with probability prec(a;—-Z^*(a;)) where Z* is the location 
of the cluster-head nearest to x. 

In order to formalize the problem of the economic opti- 
mization of the proportion of the two types of devices, let 
us assume a regular repartition of cluster heads on the plane. 
A simple model consists in taking them to be repartitioned 
on a regular, say, triangular grid with some density Ac. This 
means that Ac = 4/(i^V3), where L is the distance between 
two adjacent cluster-heads. Note that maximal di stance to a 
nearest cluster head is equal to i?max(Ac) ~ 4/ \/ Ac3\/3. As 
as in the previous section we model the traffic of packets sent 
be the sensors to the cluster-heads (who act independently) by 
Poisson rain model of events that is assumed to be stationary 
both in time and on the whole plane M^. To further simplify 
the model, we assume that at each point x in space at least one 
cluster-head has to achieve p{x) > D. This is an upper bound 

^In other words, to maximize the total capacity it is optimal to receive only 
packets whose received power is larger than some threshold. 



on Ac; in reality, a packete that is lost by a cluster-head, may 
still be captured by another cluster-head. However, this upper- 
bound is sufficient to numerically demonstrate large savings 
of the hybrid approach. 

Consider the following problem: minimize the cost of the 
network by unit area C — XgCs + AcCc given some minimal 
intensity of received information AgAc + /o(-Rmax) > D, 
where Ac Ac is the density of information captured directly 
by the cluster-heads and /o(-Rmax) is the lowest density of 
information that can by obtained from the sensors given max- 
min (maximizing coverage) admission policy. 

In order to solve this problem, given D and Ag, A^, we take 
the max-min poHcy d^^xm < I3.2> with D{x) = l(|a;| < R) 
and find the maximal radius Rmaxm^ such that the constant 
P-maxm Obtained by the policy D{x) = l{\x\ < R„,axm) 
on i?(0, i?„Q^„j) is equal to D ~ AgAc (cf. Example \3.5i . 
Note by Proposition 13.41 that for all x, p{x) > p{x) > 

PmaxnS^rnaxm) = D - XeXc- This meanS that taking Sinax = 

Rmaxm'^ i-C-, Ac = 4/ (i?^„^„3\/3) is Sufficient for AcAc + 
/5(-Rmax(Ac)) > D. Having calculated Ac = Ac(As) we express 
the cost of the network C = C{Xc/Xs) as the function 
the proportion between the intensity of cluster-heads and the 
sensors. Finally, we look for its maximal value. 

IV. Numerical results 

We will give now some numerical examples. We consider 
the canonical traffic scenario described in Section III-AI with 
SINR threshold 7=1, path-loss model L{x) = L{\x\) = 
with K = 10"5•^ ri = 3.3. 

a) Single cluster-head scenario: In this part we consider 
a scenario with a single cluster-head at 0. Figure □(left) shows 
the the quality of approximations given in Proposition 13.31 

Next we compare maximum radii of different admission 
policies. In addition to d^^axm d*, defined in Sec- 
tion IIII-BI and Section IIII-CI respectively, we introduce the 
coverage-optimal deterministic (COD) policy. It is defined as 
dcooix) — l(|a;| < Rcod) where Rcod is the maximum 
radius such that p[Rcod) > D under poHcy dcoD- 

We can see from the results in Figure |4] (left) that Rmaxm ~ 
Rcod, hence that a deterministic policy with a well-chosen 
radius provides almost as good coverage as the max-min fair 
policy. In addition, we see from Figure El (right) that COD 
policy is more efficient. A more detailed discussion on a 
tradeoff between efficiency and fairness is out of the scope. 

Figure |3] (right) shows the total intensity U of information 
obtained when the admission policy accepts all the packet 
within a given radius. Optimal policy radius (maximizing 
the total throughput) can be deduced from this plot. Two 
marked radii correspond to policies , d .We see that the 
policies d* , d well approximate the optimal. We also see on 
Figure H (right) that maximizing capacity requires admission 
region R* that is much larger than R„iaxm ^^'^ Rcod, having 
significantly smaller p(x) < p (x). 

b) Economic optimization: We now look at the economic 
aspects, described in Section UlI-DI Figure |5] (left) shows the 
required density of cluster-heads in the hybrid network in 
function of the density of transmit-only sensors, given the 
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Fig. 3. Exact value and approximations of Prec (left)- Maximization of the 
total intensity of information (right). 
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Fig. 4. Radii of different admission policies (left). The exact profiles of 
p{x) for different admission policies (right) when applying different admission 
policies explicitly found using the lower-bound p (this explains why the max- 
min policy does not give exactly p(x) = const for x within the admission 
distance). 




Cost CH / Cost SN 



Fig. 5. Constitution of the hybrid network (cluster-heads CH and transmit- 
only sensors SN) with minimal cost, assuring min^ p(x) = 0.75 (right). 
Network cost gain in function of the device cost per unit ratio (right). 

minimal value min^^ p{x) = p(i?max) = 0.75. Figure |5] (right) 
shows the network cost economization in function of the 
cluster-head/sensor unit cost ratio. We can see that even when 
a price of a cluster-head is only slightly higher than a price of 
a transmit-only sensor, we can achive significant savings. On 
a contrary, using the naive policy, cost savings are negligable. 

V. Implementation issues and concluding remarks 

In this paper, we analyzed hybrid sensor networks consiting 
of transceiving and transmit-only sensors. We presented a 
detailed mathematical model of a physical and MAC layer of 
the network. Using this model we derived the optimal packet 
admission policies for cluster-head that maximize coverage or 
total throughput. Also, using the model we demonstrated how 
much the dollar-cost of a sensor network can be decreased 
while maintaining the same network coverage. The MAC 
model that we developed can also be used for any non-slotted 
wireless communication channel. 



In this work we did not discuss implementation details of 
the optimal policies. However, we note that these policies can 
easily be implemented, based only on the knowledge of packet 
received power (no need to know channel attenuation function, 
sensor positions, etc). One implementation is discussed in [1], 
and can easily be generalized to the policies proposed in this 
paper. We leave details for future work. 

Also, in this work we optimized only packet admission 
policy and not the sensor density As(x), transmitting power 
P{x), coding 7(0;) nor the transmission frequency Ae(x), 
which may depend on particular sensors. However, we showed 
that even constrainted on the optimization of the admission 
policy only, we can achieve large savings and maintaint 
architecture simple. Optimizing other parameters is left for 
the future work. 
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Appendix: Mathematical modeling 

In this section we present mathematical models that are 
used to analyze the sensor network described in Section |ll] 
In particular, 

A. An Erlang 's M/D/1/1 loss system with interference 

Assume a time homogeneous, independently marked Pois- 
son point process $ = {(T„, (P„, i?„))};f^_^, where r„ 
are customer (packet) arrival epochs and (P„,iJ„) are inde- 
pendent, identically distributed (i.i.d.) marks, where P„ > 0, 
Hn > can be interpreted as, respectively, the average (over 
fading effects) power with which the nth packet arrives at 
the receiver and the actual fading state of its channel. (The 
randomness of {Pn} reflects different locations of transmitters 
and powers with which they emit packets, while the random- 
ness of {Hn} reflects the temporal variation of the channel 
conditions given fixed location of the transmitter and emitted 
power) Lets denote by A (0 < A < 00) the intensity of $; 
i.e., T„ are i.i.d. exponential r.v. with parameter A. 

We consider the following modification of the Erlang's loss 
policy. Suppose that each arrival (i.e., packet) is admitted by 
the single server of the system (i.e., starts being received by the 
receiver) if this latter is idle at the packet's arrival epoch and 
rejected otherwise. Admitted packets are received during their 
duration time B. However, packets that are rejected by the 
receiver interfere during their emission time with the packets 
that are being received. Inspired by inequality ( 12. 1> . with 
HnPn — \h\'^Prec, we will Say that the nth packet, given it is 
admitted by the receiver, is correctly received if the following 
inequality holds 

PnH„ 



> 



7, 



(A.l) 



W + 1/B J^y" I{t) dt - PnHn 
where W is some nonnegative r.v. independent of <&, 7 > 
is some constant, I{t) is the value of the following temporal 
shot-nose process at time t 
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describing the total power received at time t from all packets 
that are being sent (including the power of the packet that 
is being received; this is why we substract PnHn from I{t) 
in dAHl . 

Our goal is to calculate the frequency of 
the correct reception of packets; i.e., tt = 
limjv_+oo #{packets received among packets nr 1, . . . , /N, 
where # denotes the number. Denote by X{t) the indicator 
of the event that the receiver is busy at time t and 
X{—0) = limtioX{t). Denote by Sn the indicator of the 
event that the inequality ( lA.H holds. Let P° denote the Palm 
probability given there is a customer arrival at time 0, and 
let E'^ denote the corresponding expectation. By Slivnyak's 
theorem, under arrivals form the original stationary 
Poisson point process with an extra arrival {Tq,{Pq,Hq)) 
added at time Tq = whose mark is independent and 
originally distributed. Denote by P, E, respectively, the 
stationary probability of $ and its corresponding expectation. 
The following results is a consequence of the ergodic theorem 
(see e.g. [11, Theorem 1.6.1]). 

Proposition A.l: The limit defining tt exists P*^ almost 
surely and tt = E"[(l - X{~0))So]. 

In order to calculate tt we will first characterize under P° the 
distribution of the shot noise process I'{t) = I{t) — PqHq 
for t € [0, B) given the server is idle just before time (i.e., 
X(— 0) = 0). It will be given in terms of the conditional 
jointLaplace transform . . . ,^„) = 

E°[e-S?=i ^'(*')«' l^(-O) = 0] evaluated for Ci > and any 
finite collection of time instants ti, . . . , t„ e [0, S), n > 1. 

Proposition A.2: Consider the Erlang's M/D/1/1 loss sys- 
tem with interference. Then 

P°{X{-0) = 0} = ^ 



1 



and for ti, . . . ,tn G [0, B), n > 1, ^i, 



XB 



(A.3) 



'(ti)„ 



/'(t„)|Ji:(-o)=o(^i 

oo 
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(A.4) 
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^ (l-l(-t.o](s)) 



l-E[e-^"=i«'^'^-i[''^«)(*'^^)jds 
Note by the form of the above Laplace transform, that 
under P° and given the server is idle just before the arrival 
of the customer at 0, the shot-noise process of interference 
{I'{t) : < t < B} is driven by a non-homogeneous, double- 
stochastic Poisson process with intensity equal to A on the sum 
of the intervals {-B, -T] U (0, B) and elsewhere, where T 
is exponential random variable with parameter A. 

Proof: Note that ( IA.3t follows directly from the Erlang's 
loss formula (see e.g. [12, equation (81), p. 71]). In conse- 
quence, A/(l + XB) is the intensity of the point process of 
arrivals of packets that are accepted by the receiver In order to 
prove the remaining part of the proposition, lets define by T{t) 
the time that elapsed form the last moment before t when the 
receiver was busy; i.e., T{t) —t — sup{s : s < t, X{s) = 1}. 
We will first show that under P° and given X{—0) = 0, the 
variable T — T(0) is exponential with parameter A. Indeed, 
for M > 0, by the Neveu exchange formula (see e.g. [11, 



dt 



Section 1.3.2]), 

P°{r(0) > u,X{- 



0)^0} = P°{T(0)>ii} 



1 



1 + AB 



l{T{t) > u)$(dt) 



(A.5) 



where E° corresponds to the Palm probability, given a packet 
arriving at time is accepted by the receiver, is the next 
arrival time after of a packet accepted by the receiver, and 
the integral with respect to $(dt) denotes the sum over all 
arrival times of the process $. It is easy to see that under 
E° we have T{t) = for i e (0, B] and T{t) ^ t - B for 
t G (i?,rf]. In the interval {B.T^] point process $ has just 
one point, namely , and thus the integral in \A.5\ reduces 
to l(T(rf ) > It) = l(Tf - B >u). The distribution of the 
points of $ in (0, cx)) is not influenced by the condition that 
the server is idle just before and that there was an arrival 
at 0. Thus, under Ej^, as well as under E, it is equal to the 
distribution of points of the the original Poisson point process. 
Thus, due to the lack of memory of the exponential inter- 
arrival r.v., the variable Tf — i? is exponential with parameter 
A, which completes this part of the proof. 

Note now, that the packets which contribute to /'(i) = 
I{t) - PqHo for t G [0, B), given T(0) = T and Tq = 0, 
arrive only during the time intervals {—B, — r)U (0, B). (Note 
that P°( there is an arrival at -T) = if T > 0.) Note also, 
that this set is disjoint with the set 5 = [T*, -B - T(0)] U 
[-T(0),0], where T* = sup{r„ : T„ < T(0) - B,T„+i < 
T{0)~B, Tn+i—Tn > B}. This latter set is a random stopping 
set (i.e., for a given realization of the point process the set 
5(0) is invariant with respect to any modification of the points 
of in K\5(0); see e.g. [13]). Also, X{-0) and r(0) depend 
only on the configuration of points of (p in 5(0). Thus, by the 
strong Markov property of the Poisson point process (see [13]), 
given X{-0) = 0,T = T(0), the distribution of amvals in 
{—B, — r)U(0, B) is equal to the original distribution of points 
of the independently marked Poisson point process $ taken on 
this region, and hence iAA\ holds. This completes the proof. 

■ 

Before we give an explicit formula for the frequency tt 
of the correct reception of packets for our loss system with 
interference in the case of Rayleigh fading, we will calculate 
the conditional Laplace transform of the integrated shot-noise 
I'b = 1/Bj^r{t)dt given P° and X(-0) = 0; i.e., 
^?^|X(-o)=o(0 = EO[e-«^^^ \Xi-0) 0] for ^ > 0. 

Proposition A.3: Suppose {Hn} are exponentially dis- 
tributed with mean 1 and {Pn} are independent of {Hn}- 
Then /:?m;,(_o)=o(0 = (0^2(0, where 



CiiO = exp(^-Ai?(l-E[^log(l + eP) 
£2(0 = exp(-AB) 



1 + XB exp ABE 



(A.6) 
(A.7) 
dt 



Proof: Note first, that the integrated shot noise 
is also a shot noise type variable. Indeed, I'g = 

^/Bj,''Ei^oHim[o,B){t-T,) - E^^oH^P^ViT,), whcrc 
V{t) is equal to B — \t\ for \t\ < B and elsewhere. 
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By the Proposition IA.2I it suffices to calculate the Laplace 
transform of I'g driven by independently marked double- 
stochastic Poisson process with intensity A on the sum of the 
intervals {—B, —T] U (0, B) where T is exponential random 
variable with parameter A. Denote 

(A.8) 



7^2 



Y,H^P^V{T,)\_B,-T]{T^)■ (A.9) 

Applying the general formula of the Laplace transform of the 
independently marked Poisson point process (see e.g. [14])) 
we get 



E*'[e-«-f«|X(-0) = 0] =cxp -A/ (l-E 



,-(,HPV(t) 



At 



(A. 10) 

where H, P are independent, generic r.v's for {i?„} and {Pn}- 
Integrating with respect to t and evaluating expectation with 
respect to the exponential rv. H we obtain E[e^?-^B] = Ci{^). 
In order to calculate ¥P[e-^^B\x[-Q) = 0], we condition 
on T = r(0) and we use similar arguments, with integral 

fB f — T^ 

Jp in (IA.10> replaced be (note that this integral is 

null if T > B). Similarly as for /jj and then by integra- 
tion with respect to the law of exponential T we obtain 
E°[e-«-fB|X(-0) = 0] = £2(0- Obviously I'b = ih + 
and variables /Jj , are independent, because they are driven 
by disjoint regions of the underlying Poisson point process. 
Thus (^) — £i(^)£2(0' which completes the proof. ■ 

Now we are able to give the main result of this Section — 
an Erlang's type formula. 

Proposition A.4: Consider the Erlang's M/D/1/1 loss sys- 
tem with interference. Suppose that {Hn] are exponential rv's 
independent of {Pn} and lets denote by Cw{0 'he Laplace 
transform of W . Then 

where the expectation is taken with respect to the random 
variable Po- 

Proof: By Proposition lA. II tt is equal to 

E"[(1-X(-0))M 
= V°{ X(-0) = } E"[,5o|X(-0) = 0] 

= P"{ X{-{)) } P"{ i/o > lW{W + /^)/Po|^(-0) - } 

= P0{ X{-0) = } E" [exp(-7(W^ + /^)/Po) |^(-0) = o" 

because iJo is exponential with mean 1 and independent 
of everything else. Conditioning on Pq, noting that W is 
independent of I'g and using Propositions IA.2I and IA.3I we 
obtain the result. ■ 
Lets introduce now to the Erlang's loss model an additional 
(external) stationary, ergodic process J{t) of interference, 
independent of W and (For example, one can think of J(t) 
as of the interference created by emitters transmitting packets 
that are not supposed to be received by our receiver due to 
some random, independent admission policy; cf. Section lA-Bt . 
Lets say that the n th packet of given it is admitted by the 
receiver, is correctly received if to the inequality (lA.U holds 



with the term \/B /^"^ J{t)At added to the denominator. 
Denote by 5'^ the indicator of this event, and it' — E"[(l — 
X{-Q))5'f^]. Denote by Jb = J{t)At and its Laplace 
transform by Cjg{£,) — E[e^^'-'^]. We have the following 
straightforward extension of the Proposition lA. 41 

Corollary A.5: Consider the Erlang's M/D/1/1 loss system 
with the external interference J{t). Under the same assump- 
tions as in Proposition IA.4I the frequency of the correct 
reception of packets exists P° almost surely and is equal to 
, E [Cw {l/Po)Ci h/Po)C2 {l/Po)Cj, (7/Po)] 

TT = . 

l + XB 

B. Sensors on the plane 

In this section we assume that the packets are emitted from 
different locations of the plane and, assuming some form 
of the attenuation function, we will obtain a particular form 
of the distribution of the powers {Pn} received at the origin, 
where the receiver is supposed to be located. (Note that this 
distribution was not specified in the previous section.) We will 
also assume some packet admission policy. 

Attenuation function: Suppose that the signal transmitted 
with some power P at the location x is attenuated on the path 
to the receiver located at (on average over fading effects) by 
the factor L{x) > 0.; i.e., the average received power is equal 
to PL{x). 

Spatial policy of packet admissions: Suppose that packets 
are emitted from different locations of the plane M^. Suppose 
moreover, that the receiver located at the origin adopts the 
following spatial admission policy. Depending on emission 
location x, it accepts the packet, independently of everything 
else, with probability d{x) (and starts receiving it, provided it 
is idle), where < d{x) < 1 is a given function of location x. 

1) Poisson rain model of events: Consider a spatio- 
temporal Poisson process {(X„,r„)} where X„ € B C 
K^, r„ e M, with intensity measure As{dx) x Xedt. The 
coordinates of the point (X„,T„) denote, respectively, the 
location of a packet emission and the time it starts. (One can 
think of emitters being born at locations Xn and time r„ just 
to emit one packet at this moment; after the transmission of 
this packet the emitter disappears.) We assume that the points 
{{Xn,Tn)} are independently marked by i.i.d. exponential 
(with mean 1) random variables Hn modeling the fading 
conditions during the transmission n. Moreover, assuming 
some admission policy d{-) we suppose that the points are 
further marked by i.i.d Bernoulli variables [/„ describing 
the admission status of the packets; i.e., P{ [/„ — 1 } — 
1 — P{ Un — 0} — d{Xn), where [/„ = 1 marks an admissible 
packet. 

We call the marked Poisson point process 4' = = 
|((X„, T„), (_ff„, [/„)) } , the Poisson rain of events with a 
given spatial admission policy d{-). We consider ^E* as the 
input to the Erlang's loss system with interference described in 
Sections [A- Al Specifically, we define P„ = P L{Xn) and take 
the admissible packet transmissions $ = {{Tn, {Pn, Hn)) : 
Un — l}n as the input to the system, whereas the total received 
power from non-admissible packet transmissions 

J{t) = UnHnPL{Xn)l[0,B){t - Tn) , (A.14) 
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A(0 
MO 



= exp -AB + AeS 



d{x) 



exp(-AB) + Xsj^ exp(^AeB ^ log(l + (,PL{x)t) A,(dx) ) di ) , 

exp (^-2KB j{l - d{x)) (l - log(l + iPL{x)) ] A,(dx) 



(A.ll) 
(A. 12) 
(A.13) 



as the external interference. Denote by 



A = A, 



d{x) As (da;) 



(A.15) 



the (temporal) arrival intensity of Poisson process of the 
packets admissible according to the spatial policy d{-). The 
following consequence of Corollary IA.5I gives the Erlang's 
type formula for the Poisson rain model. 

Corollary A.6: Consider the Erlang's loss system 
M/D/1/1 driven by the Poisson rain of events on some 
domain B with spatial admission policy d{-). Assume that 
A given by ( IA.15> is finite. Then, the fraction it' = tt'{xo) 
of admissible packets correctly received from a location 
xq E M^, given there is an emitter located there, is given by 
Corollarv IA.5I with constant Pq = PL[xa), and Ci^ C2, Cj^ 

[3, (IA.13>. where A is 



given, respectively, by (lA.lH . 
given by (IA.15> and the integrals are taken over D. 

Proof: Note that the distribution of the received 
power is equal to P{ P < a} — J d{x)l[PL{x) < 
a) As(dx)// d{y) As{dy), which is correctly defined since 
we assume \ < 00. Then, formulas ( lA.lH . (IA.12> follow, 
respectively, from iA.6\ . (IA.7> . Next, note that $ and J{t) 
are independent; this is a consequence of the independent 
thinning of the Poisson process of all packet emissions. 
Moreover, the integrated interference Jb, given by the formula 
Jb = - Un))HnPL{Xn)V{Tn) IS a shot-noisc type 

random variable (cf. the proof of Propositon lA.3> . Its Laplace 
transform Cj^{£,) is known explicitly and given by ( IA.13> . ■ 
2) Fixed arbitrary locations of emitters: Suppose emitters 
are fixed and located at {Xi}. This case can be seen as a 
special case of the Poisson-rain of packets, with purely atomic 
spatial density measure As{D) = if{Xi : Xi e D}. Then the 
integrals J{. . 



, ) dx in formulas ( IA.ll> -( lA.13> take form of the 
respective sums (■• ■)• Moreover, if the spatial repartition 
of Xi is sufficiently dense, then these atomic measures can 
be reasonably "smoothed" leading to approximative integral 
formulas. In particular, if the repartition is dense and uniform 
(in empirical sense) then the sums can be approximated by 
integrals with respect to the Lebesgue's measure A(da;) = 
As da; with A, = #{X,}/|©| where |e| is the surface of D. 

3) Bounds: In this section we will give some simple bounds 
for the frequencies of successful receptions of packets. Denote 
7. = 7/{PLix)). 

Corollary A.7: Under the assumptions of Corollarv lA. 61 we 
have 



< 7r'(a;o) < 



XoJ 



Proof: Note first that £2(0 > ^liO- This can be 
verified directly comparing formulas (IA.6> and (IA.7> . but a 
simple probabilistic argument can be used as well; remind 
that that £i{£,) is the Laplace transform of /]j given by ( IA.8> . 
whereas £2(0 is the Laplace transform of Ig given by (IA.8> . 
Than, the lower bound of Corollary IA.7I follows immediately 
from 



(lA.lH and (IA.13> . In order to get the upper bound, it is 
enough to observe that £2 < 1 and to take Cj^ with factor 2 
in the exponent of ( IA.13> replaced by 1 . ■ 
Note that the upper bound in Corollary IA.7I consists in taking 
no interfering arrivals before reception of a given packet, 
whereas the lower bound consists in assuming the uncondi- 
tional Poisson process of such arrivals (T = in the proof of 
Proposition IA.2l l. 

References 

[1] B. Radunovic, H. L. Truong, and M. Weisenhorn, "Receiver architectures 

for transmit-only, UWB-based sensor networks," 2005. 
[2] A. Akyildiz, "ESRT: Event-to-sink reliable transport in wireless sensor 

networks," in Pmc. of ACM Mobihoc, 2003. 
[3] B. Otis, Y. Chee, and J. Rabaey, "A 400uW-RX, 1.6mW-TX super- 
regenerative transceiver for wireless sensor networks," in Proc. of IEEE 

International SoUd-State Circuit Conference, 2005. 
[4] I. O'Donnell, "UWB transceiver project update," in Proc. of BWRC 

Summer Retreat, June 2003, poster Session. 
[5] S. Sheng, L. Lynn, J. Peroulas, K. Stone, and I. O'Donnell, "A lowpower 

CMOS chipset for spread-spectrum communications," in Proc. of IEEE 

International Solid-State Circuit Conference, 1996, pp. 346-347. 
[6] B. Liu and D. Towsley, "On the coverage and detectability of wireless 

sensor networks," in Proc. of WiOpt, 2003. 
[7] O. Dousse, R Mannersalo, and P. Thiran, "Latency of wireless sensor 

networks with uncoordinated power saving mechanisms," in Proc. of 

ACE Mobihoc, 2004, pp. 109-120. 
[8] S. Janson, "Random coverings in several dimensions," Acta Math., vol. 

156, pp. 83-118, 1986. 
[9] H. Koskinen, "On the coverage of a random sensor network in a bounded 

domain," in Proc. of 16 ITC Sepecialist Seminar, Antwerp, Belgium, 

September 2004, pp. 11-18. 
[10] D. Tse and P. Viswanath, Foundamentals of Wireless Communication. 

Cambridge University Press, 2005. 
[11] F. Baccelli and P. Bremaud, Elements of Queueing Theory, ser Serie: 

Applications of Mathematics. Springer Verlag, 2002, second edition. 
[12] R. Wolff, Stochastic Modeling and the Theory of Queues. Prentice-Hall, 

1989. 

[13] S. Zuyev, Strong Markov Property of Poisson Processes and Slivnyak 
Formula, ser. Lecture Notes in Statistics. Springer, 2006, vol. 185, pp. 
77-84. 

[14] D. Stoyan, W. Kendall, and J. Mecke, Stochastic Geometry and its 
Applications. Chichester: John Wiley and Sons, 1995. 



1 + AB - ■ - l + \B 

where C and A are given by and (IA.15> . respectively. 
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